Spin transport in the Neel and collinear antiferromagnetic phase of the two 
dimensional spatial and spin anisotropic Heisenberg model on a square lattice 
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We analyze and compare the effect of spatial and spin anisotropy on spin conductivity in a 
two dimensional S=l/2 Heisenberg quantum magnet on a square lattice. We explore the model 
in both the Neel antiferromagnetic (AF) phase and the collinear antiferromagnetic (CAF) phase. 
We find that in contrast to the effects of spin anisotropy (A) in the Heisenberg model, spatial 
anisotropy (rf) in the AF phase does not suppress the zero temperature regular part of the spin 
conductivity, a reB (uS), in the zero frequency limit - rather it enhances it. We also explore the finite 
temperature (T) effects on the Drude weight, Ds(?7, A, T), in the AF phase for various spatial and 
spin anisotropy parameters. We find that Ds goes to zero as the temperature approaches zero. At 
finite temperatures (within the collision less approximation) enhancing spatial anisotropy increases 
the Drude weight value and increasing spin anisotropy decreases the Drude weight value. In the 
CAF phase (within the non-interacting approximation) the zero frequency spin conductivity has a 
finite value for non-zero values of the spatial anisotropy parameter. In the CAF phase increasing 
the spatial anisotropy parameter suppresses the a reg (cj) response at zero frequency. Furthermore, 
we find that the CAF phase displays a spike in the spin conductivity not seen in the AF phase. 
Inclusion of the smallest amount of spin anisotropy causes <7 re9 (aj) to develop a gap in the spin 
conductivity response of both the AF and CAF phase. Based on these studies we conclude that 
materials with spatial anisotropy are better spin conductors than those with spin anisotropy both at 
zero and finite temperatures. We utilize exchange parameter ratios for real material systems such as 
SrZnVO(P04)2 (spatially anisotropic) and La2Ni04 (spin isotropic) as inputs to the computation 
of spin conductivity. 

PACS numbers: 75.10.Jm, 75.30. Ds, 75.40. Gb, 75.76.+J 



I. INTRODUCTION 



The concept of spin transport is both intriguing and 
exciting [l|, Q . The importance of studying spin trans- 
port is related to the field of spintronics which involves 
the study of active control and manipulation of spin de- 
grees of freedom in a solid state material [j-Q. The 
fundamental quantity of interest is spin conductivity as- 
sociated with the flow of spin current. Charge currents 
in two-dimensional high-mobility electron systems with 
Rashba spin-orbit coupling can generate spin currents 
However, power dissipation in spin systems in 
which spin transport is accompanied by charge transport 
limits the possible application to state-of-the-art spin- 
tronic devices. In nonitinerant quantum systems the dis- 
sipation problem is reduced. This is a major technologi- 
cal motivation behind exploring the fundamental physical 
processes of spin transport in insulating magnets [lOj. 

The definition of spin current is a controversial topic 
[Tl] - [l^ | . The issue arises in the context of spin transport 
in semiconductors where spin is not conserved because 
of the presence of spin-orbit interaction. For insulating 
magnets the spin current operator definition, Ji^j, has 
been investigated by Schiitz et al. It has been shown that 
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for a collinear spin configuration (the case considered in 
this paper) magnetization transport is appropriately de- 
scribed by the scalar current density operator defined in 
Eq. [T5] Therefore issues associated with the conceptual 
definition of spin current are of no concern for the spin 
conductivity analysis carried out in this paper. Possible 
experimental set up to measure spin conductivity have 
been proposed in Ref. UH. 

Spin tr ansp ort phenomena has been explored in ferro- 
magnetic [la, [I?} and antiferromagnetic (AF) insulators 
[18h24| . Low-dimensional magnetic systems have been 
investigated using the spin anisotropic Heisenberg AF 
model for both S=l/2 and S=l. The 2D & 3D Heisen- 
berg model [25| and 2D XY magnets [26| have also 
been analyzed for spin transport. Additional work has 
been done on the effects of Dzyaloshinskii-Moriya inter- 
actions on spi n Heisenberg magnets in one- and two- di- 
mensions [27[. Present analysis on the models explored 
till date show unconventional ballistic spin transport at 
finite temperatures. For the two-dimensional (2D) AF 
the regular part of the spin conductivity remains finite in 
the dc limit (for the isotropic Heisenberg system) [25j |. 
at zero temperature, while for the ID case the regular 
part of the conductivity is suppressed at low frequen- 
cies [ill, [2(|. Two-magnon processes contribute to the 
magnetization transport and can be attributed to the mo- 
tion of quasi-particles as magnons (2D) or spinons (ID). 

In this paper we study the effects of spin con- 
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FIG. 1. (Color online) The spatial and spin anisotropic 
Heisenberg model considered in this paper (see Eq. [TJ. Ex- 
change interactions along the x- and y- direction are given by 
J x and J y respectively. The z- component spin anisotropy is 
given by A. The above image is for the model in the anti- 
ferromagnetic phase. The model is studied both in the Neel 
antiferromagnetic and collinear antiferromagnetic phase. The 
external magnetic field B z (x,t) is directed along the z- direc- 
tion and is spatially inhomogeneous along the x- direction. 
The inhomogeneous magnetic field drives a spin current, j x , 
along the x- direction. In this article we explore the effects of 
spatial and spin anisotropy on the spin current ] x . 



ductivity on a spatial and spin anisotropic S = 1/2 
Heisenberg model (see Fig. [T]). The model is not 
purely of academic interest since there are several 
vanadium phosphate material systems (Pb2VO(P04)2, 
SrZnVO(P0 4 ) 2 , BaZnVO(P0 4 ) 2! and BaCdVO(P0 4 ) 2 ) 
in which extensive bandstructure calculations show a spa- 
tially anisotropic exchange interaction along the x- and 
y- direction [28|. For isotropic exchange interactions we 
consider La2NiC>4 as the model system [2^. We utilize 
these real material exchange parameter ratios as inputs 
to compute the spin conductivity response in the Neel 
AF phase. We discuss the temperature (T) dependence 
of the spin Drude weight, T>s(r), A, T), on spatial and 
spin anisotropy within the collision less approximation 
scheme for the magnons. We also analyze the spin con- 
ductivity behavior of the previously unexplored collinear 
antiferromagnetic (CAF) phase. 

The main results of this paper can be summarized as 
follows (i) In contrast to the spin anisotropic Heisenberg 
model, spatial anisotropy in the AF phase does not sup- 
press the regular part of the spin conductivity in the zero 
frequency limit. Rather it helps to enhance the spin con- 
ductivity. Compared to spatial anisotropy we find that 
spin anisotropy strongly suppresses the zero frequency 
weight of the regular part of the spin conductivity (see 
Figs. El HI and [S|), (ii) as expected the finite temperature 
Drude weight in the AF phase vanishes in the zero tem- 
perature limit (this is the same behavior displayed in the 
spin anisotropic 2D Heisenberg model). Increasing spa- 
tial anisotropy increases the Drude weight value. How- 
ever, increasing the spin anisotropy parameter decreases 



the Drude weight value (see Fig. [6J. The Drude weight 
results are obtained within the non-interacting (collision 
less) magnon approximation (see Section IIVBI for dis- 
cussion) and (iii) spin conductivity in the CAF phase 
has a non vanishing zero frequency contribution. In con- 
trast to the AF phase, spatial anisotropy in the CAF 
phase decreases the zero frequency regular spin conduc- 
tivity value. The CAF phase also displays a spike in the 
spin conductivity. Based on the above analysis we con- 
clude that materials with spatial anisotropy are better 
spin conductors than those with spin anisotropy in both 
the AF and CAF phase at zero and finite temperatures. 

This paper is organized as follows. In Section [H] we 
introduce the model Hamiltonian for the AF phase and 
derive the Fourier transformed Dyson-Maleev (DM) ver- 
sion of the Hamiltonian upto quartic interactions. The 
vertices for the quartic interaction are listed in Ap- 
pendix [AJ In Section IIIII we state and derive the ex- 
pression for the spin current operator for our model. 
In Section IIVI we utilize the many body Green func- 
tion formalism to compute the spin conductivity. In 
Section IIV Al we present and discuss for the AF phase 
the non-interacting (Section IIV A 1[) and interacting spin 
conductivity (Section TlV A 2[) results. In Section llVBI we 
discuss the anisotropy and temperature dependence re- 
sults of the drude weight contribution. In Section IIV CI 
we introduce the model Hamiltonian for the CAF phase 
and present the results of the spin conductivity within 
the non- interacting approximation. Finally, in Section IVl 
we summarize the main findings of our work. 



II. MODEL HAMILTONIAN 

The Hamiltonian for the 2D Heisenberg model on the 
square lattice with both directional anisotropy of ex- 
change couplings and spin anisotropy is defined as 

H = £ U\{S+S- + ^ + S~S+ +X ) + ASfSf +x ] 
+ E Ul(S+Sr +y + S r S + +y ) + AS?Sf +y ], (1) 

The exchange interaction in the x-dircction is denoted 
by 3 X and along the y-direction by 3 y . The spin vector, 
Si , is defined on every lattice site where the index i runs 
over all lattice sites and i+x and i +y runs over nearest 
neighbors in the x- and y- direction respectively. S^ and 
S~ represent the raising and lowering spin operators at 
site i. Sf represents the z-th component of the spin vector 
Si . The spin anisotropy parameter is given by A and the 
spatial anisotropy parameter by ?y = ^j-. The magnetic 
model considered here supports three types of classical 
ground-state configurations - the Neel (7r,7r) AF state and 
columnar (ir,0) or row (0,7r) CAF phase. The AF state 
survives for r\ > & A > 1. The CAF phase for 77 < & 
A > 1 . In Table Q] we list examples of material systems 
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with exchange interaction ratios and spin anisotropies of 
interest in this paper. 

We study the above Hamiltonian using the DM trans- 
formation. In the DM representation the spin operators 
are replaced by the bosonic operators as 

S i * = 5-o i t o i ,S , r = v^So i t , 

St = V2S^l-^)a i , (2) 

for the up-spin (A) sublattice and by 

S*~ = -S + b\b,, Sj = V2Sbj, 

S+ = yffib] (l - b -§) , (3) 

for the down-spin (B) sublattice. aj(a^) and bj(b,) rep- 
resent the creation (annihiliation) operators on site i in 
the A or B sublattice respectively. The Fourier transfor- 
mations for each operator and b, is defined below 

* = y|E e " ik ' Ri ^^ - y|E« +ik - R ^- ( 4 ) 

k k 

where N is the number of lattice sites and the momen- 
tum k is restricted to the magnetic Brillouin zone. Uti- 
lizing the definition of the Fourier transformations, the 
DM transformation in Eqs.[2]and[31 and the Bogoliubov 
transformations 

a k = u k Q k + u k /3 k , fo k = WkA< + vical, (5) 

we can diagonalizc the quadratic part of the Hamiltonian. 
The original Hamiltonian, Eq. [TJ can then be written in 

I 



k-space as 

H DM = E + H x + : V DM : (6) 

where Hdm is the DM transformed Hamiltonian. The 
classical energy (Eo) and the linear spin wave theory 
(LSWT) hamiltonian (Hi) is given by 

E = -NAJ x (l+ V )S 2 a 2 (S), (7) 
Hi = ^ fiflk(<4a k + (8) 

k 

The coefficients Uk and v k are given by, 

Mk = \r5? ■ Wk = " S5?l(7k) \/^f ■ (9) 

The LSWT dispersion (ilk) and the defintion of 7k is 

Ml k = 2AJ X (1 + ?? )5a(5)e k , (10) 

e k = ^/l-7S/A 2 , (11) 
7k = — - — -(cosk x +1] cos k y ). (12) 

0- + V) 

The Oguchi correction [3(| cx(S) is given by 

a(5) = 1 + ^ ;r = 1 -|E £ k- (13) 

k 

The normal ordered DM quartic interactions, 'N dm'-, can 
be expressed in terms of the a and /3 bosons as 



: V DM ■ = N Sg ^ + 2 - 3 - 4)(U (1) a+a+a 3 a 4 + U (2) a+/3 2 a 3 a 4 + V® atf <4 ft a* + U (4) a+a 3 /3+£ 2 

(1234) 

+ F (5) /3 4 + a 3 /32/3i + V^ptp+a+fa + V^a+a+(3+f3+ + 7 (8) /?i/3 2 a 3 a 4 + V {9) p+ p+ ■ (14) 



where G is a reciprocal-lattice vector and the symbols 1, 
2, 3, and 4 stand for wavevectors ki, k 2 , k 3 , and k 4 re- 
spectively. The expression for the vertex functions, V^ m \ 
with m=l,2,...,9 are stated in Appendix A. The interac- 
tion vertices include the effect of spin anisotropy through 
the A parameter and is a generalization of the vertices 
stated in Ref. HH 

III. SPIN CURRENT OPERATOR 

It has been shown that the presence of an inhomoge- 
neous magnetic field can drive a spin current (25j . We 
consider a time dependent magnetic field directed along 
the z-direction with a spatial modulation along the x- 
axis, B=B z (a;,i)z (sec Fig. [1]). The Hamiltonian, Eq. Q] 



in the presence of time dependent external magnetic field 
can be rewritten as 

H{t) = H -J2s z (x)B z (x,t). (15) 

X 

Using the continuity equation for the S z component of 
the spin we can write down the basic definition of the 
spin current operator [25[ 

j i -> j = ±J ij (S+Sr-SrS+). (16) 

Since there is no symmetry breaking field present along 
the y-dircction we do not expect any currents to flow 
in that direction. We therefore focus on the longitudinal 
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TABLE I. Exchange anisotropy parameter ratios used in the 
computation of spin conductivity. The spatial anisotropy pa- 
rameter is given by rj and the spin anisotropy parameter by 
A. See Figs. O-E] for results. 



Material 


Type 


r/=J y /J x ,A = 1 


La 2 NiO,£ 


Spatially isotropic 


1, 1 




Spin isotropic 




SrZnVO(P0 4 )^ 


Spatially anisotropic 


0.7, 1 




Spin isotropic 





a Ref. 
b Ref. 



spin transport current j x of the S z component of the mag- 
netization. Using Eq. [16] we obtain the following expres- 
sion for the spin current operator in real space (within 
LSWT) 



jxo — SJi 



J2( aibl +* ~ a t h i+x) + J2^ a i+x - hai +x ) 



leA 



leB 



(17) 

We then Fourier transform the spin current operator uti- 
lizing Eq. |4] to obtain 



3x0 = 2J x Sa(S) 2J sin(fe x ) 

k 

-2J x Sa(S)]Tsin(fc x ) 



— (a+a k + /?+/3 k ) 



tk 



.(18) 



Spin conductivity in the long wavelength limit can 
be computed using the general definition provided in 
Ref. HI We have 



2 (~ Kx) - Ax^(q = o,^) 
i{uj + i0+) 



(19) 



where A xx is the longitudinal spin current correlation 
function 



A, 



dfe^+ i0+ )(b*(q,i)^(-q,0)]>, (20) 



hN J 

and (K x ) is the spin- flip operator 



i-\-xl 



(21) 



g is the gyromagnctic ratio, /is is the Bohr magneton, 
and h is the reduced Planck constant. Based on the above 
definitions the real part of the spin conductivity can be 
written as 



$le[cT X x}=D s 5(u J )+a re 9(u;), 



(22) 



where the spin Drude weight, Ds, associated with the 
singular part of the conductivity is defined as 



D, 



(gfi B ) 2 {(-K x ) - Ke[A M (q = 0,w ->■ 0)]}, (23) 



and the regular part, a r xx g (uj) is 

3m[A m (q = Q,u))] 



(24) 



The Drude weight, Ds, is the zero- frequency contri- 
bution for the real part of the spin conductivity. In our 
model the Drude weight is dependent on T, 77, and A (sec 
Fig.[6j). A finite value of Drude weight indicates ballistic 
transport. At zero temperature Ds=0 is a signature of a 
spin insulator and Ds > of a spin conductor [ID, [33| . 
This is similar to the classification scheme for charge con- 
ductivity [Hj]. In the next section, Section [TV] we use 
the defintion of the spin current operator and the many 
body Green function formalism description to compute 
the effects of corrrclation on the spin conductivity. 
IV. SPIN CONDUCTIVITY 

In this section we derive explicit expressions for the 
regular part of the spin conductivity, a xx 9 {ui). We pro- 
ceed by defining the magnon propagators 

G aQ = -i(0|Ta k (t)4(0)|0), (25) 
G pp = -i<0|T/8£(t)Ac(0)|0>, (26) 

where T is the time ordering operator and |0) is the 
ground state wavefunction. The bare Fourier trans- 
formed propagators for the a and f3 magnons in the ab- 
sence of interactions are then given by 



G (o) 



1 



G 



(0) 



1 



-n k + *o+' ^ 



£0+' 



(27) 



where fi k is the LSWT energy dispersion. The a magnon 
propagators are represented by lines with single arrows 
and the /3 magnon propagators are represented by lines 
with double arrows in the Feynman diagram. The def- 
inition of a xx 9 {ui) can now be rewritten in terms of the 
magnon propagators as 

<J r x l 9 {u) = Sm[A„(q^0, M )] = Qm[G(q = 0,^)] 
(ffMs) 2 w w 

(28) 

To evaluate the above we define the time-ordered spin 
current correlation function 1251 



G(t) = -—<0\Tjx(t)3x(0)\0>, 



(29) 



where the symbols have their usual meaning. Now uti- 
lizing the definition of j x from Eq. [15] and following the 
steps outlined in Ref |25j we can write the propagator as 



[2J,5a(5)] 2 
hN 



E k , k ' sin t;; fc - n kk ,H, (30) 



where LI kk / (i) is the two-magnon propagator and is de- 
fined as 

n kk '(i) = -z(0|ra k (t)/3 k (t)a k ,(0)/3 k ,(0)|0). (31) 



5 



The two-magnon propagator can be re-written in terms 
of the single magnon propagators and a vertex function 

n kk < (w) = i fZo #G QQ (k, w + o/)G^(k, a/)r kk , ( W> w'). 

(32) 

The vertex function is denoted by r kk /(w,w') and satis- 
fies the Bethe-Salpeter equation given below [35| 

iw(«,«0 = - t^^Ek, rL 1^< kl] >>i) 

x G QQ (ki,o; + w')G (3/? (ki,a;i)r klk (w,wi). (33) 

In the above expression the four-point vertex is repre- 
sented by V kk k k (c*/, wi) and contains all the irreducible 
interaction terms. Integral equations can be set up to 
solve for the expression for the spin conductivity with in- 
teractions (ladder approximation) at zero temperature. 
We utilize these integral equations to numerically com- 
pute the effects of spatial and spin anisotropy in the AF 
and CAF phase. 

A. Antiferromagnetic phase 




FIG. 2. (Color online) Two-magnon scattering ladder dia- 
gram contributions. The propagator with a single arrow cor- 
responds to the a magnon. The propagator with two arrows 
correspond to j3 magnons. The magnon pairs undergo re- 
peated scattering as shown by the dashed line. The infinite 
sum of bubbles give rise to the interaction contribution in the 
spin conductivity. 



singularity is cured by introducing interactions as shown 
in the next section, Section HV A 21 

2. Spin conductivity: ladder approximation 



1. Spin conductivity: non-interacting 

We first investigate the behavior of the regular part of 
the spin conductivity in the zero frequency limit. This 
is equivalent to the replacement of the Green function 
with the non-interacting propagator and setting 



V, 



kk ilci k 



(a/,wi)=0. The magnon propagator in the non- 



interacting case is then written as 

C (Q) M _ i^SajS)} 2 v sin 2 (fc x ) 1 
k K 

(34) 

Based on the above expression we can obtain a re9 (uj) in 
our non-interacting case as 



reg 
xx 



{9 VI 



(I + rj) 2 ui N 



E 



s\VL 2 (k x ) 



8{u - 2e k ), 



(35) 
with 



where h is the Planck constant, ui 
ttmax = 2J X (1 + rj)Sa{S)/h. 

In Fig. [3] we display the results of the regular part of 
the non- interacting spin conductivity, a^(oj), in units of 
(g/is) 2 /h for the spatially anisotropic but spin isotropic 
(A = I) 2D Hcisenberg model in the AF phase. The 
plot is for 77 = J y /J x values of 1 (black circles), 0.7 (red 
stars), and 0.5 (blue diamonds). The choice of anisotropy 
parameters and the corresponding example material sys- 
tem is listed in Table HI From our calculations we find 
that decreasing rj (increasing spatial anisotropy) leads to 
an increase in the regular part of the spin conductivity. 
The finite value of the spin conductivity in the zero fre- 
quency limit is in complete contrast to the results of the 
spin anisotropy only model. In that model the smallest 
amount of spin anisotropy introduced a gap. The diver- 
gence near ui=2 is due to the Van-Hove singularity. This 



We obtain the corrections to the spin conductivity 
due to interactions within the ladder approximation of 
two-magnon scattering (see Fig. [2]). Within this ap- 
proximation scheme we replace the four-point vertex, 
V kklkik (u/, u>i) by the first order irreducible interaction 

part given by vj^ k k . The expression for this vertex is 
stated in the Appendix |A"1 Now following the procedure 
listed in Refs. [U and[H we can set up an algebraic solu- 
tion of coupled integral equations based on the following 
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FIG. 3. (Color online) Regular part of the non-interacting 
spin conductivity, a^% 9 (uj), in units of (g/is) 2 /h for the spa- 
tially anisotropic but spin isotropic (A = 1) 2D Heisenberg 
model in the antiferromagnetic phase. Here u = uj/ilmax 
where £l m ax = 2^(1 + rj)Sa(S)/h. The plot is for 77 = J y /J x 
values of 1 (black circles), 0.7 (red stars), and 0.5 (blue dia- 
monds). The dashed lines are a guide to the eye. Decreasing 
77 leads to an increase in the regular part of the spin con- 
ductivity. The divergence near £3=2 is due to the Van-Hove 
singularity. 



G 



decoupling schemes 



Ek sin (M7k-k'3k = 



_ sin(ki^) 



J2 k siu- 2 (k x )gu. 



(36) 
(37) 



Repeated applications of the above decoupling equations 
to the Bethe-Salpeter equation (Eg. [33]) and the magnon 
propagator leads to the following expression for the regu- 
lar part of the spin conductivity in the interacting ladder 
approximation 



reg 



K = 



2(l+ V )Sa(S) ■ 

To rewrite the above we utlizcd the definition below 



sin 2 ^) 



1 



(39) 



In Fig. 2] we display the results for the regular part 
of the interacting spin conductivity, a xx 9 (u;), in units of 
(g/^s) 2 /h within the ladder approximation for the spa- 
tially anisotropic but spin isotropic (A = 1) 2D Hciscn- 
berg model. The plot is for r\ = J y /J x values of 1 (black 
circles), 0.7 (red stars), and 0.5 (blue diamonds). Here 
u = u>/Cl max where Vt max = 2J X (1 + r))Sa(S)/H. As in 
the non-interacting situation decreasing -q leads to an in- 
crease in the regular part of the spin conductivity. In 
contrast to the spin anisotropic Heiscnbcrg model, the 
regular part does not vanish in the spatially anisotropic 
model in the uj — > limit. The Van-Hove divergence 
near the 2=2 point is cured after inclusion of interac- 
tions. Also decreasing 77 leads to a shift in the peak value 
towards lower frequency. 

In Fig. [5] we show the results of the regular part 
of the spin conductivity, a xx 9 (uj), in units of (g/^s) 2 /h 
within ladder approximation with both spin and spatial 
anisotropy. Here uj = ui/Q max where fl ma x = 2J X A(1 + 
r))Sa(S)/h. The plot is for rj = 0.5 & A = 1 (blue dia- 
monds), T] = 0.5 & A = 1.001 (red stars), 77 = 0.5 & A = 
1.01 (green triangles), and 77 = 1.5 & A = 1.01 (black cir- 
cles). Inclusion of the smallest amount of spin anisotropy 
suppresses the regular part of the spin conductivity and 
has a stronger effect than spatial anisotropy. Decreas- 
ing r\ from 1.5 to 0.5 but keeping the spin anisotropy 
constant (A = 1.01) simply increases the magnitude of 
a xx 9 (oj) and leads to a shift in the peak. There is no ef- 
fect on the gap above which spin conductivity is non-zero. 
The conductivity gap is controlled by the spin anisotropy 
parameter only. 

The dc conductivity is given by the zero frequency limit 



of the regular part and is defined as, a 



DC ■ 



reg 



0). A non-zero value of a DC is an indicator of a non- 
ideal spin conductor. Within our spatially anisotropic 
model in the AF phase we find several parameter ranges 
where the system is not a spin insulator rather a non-ideal 
conductor. Based on the above calculations we conclude 
that a spatially anisotropic quantum Heiscnbcrg AF is a 
better spin conductor than a spin anisotropic system at 



h (l+, ; )2A^2^' m ll+ K [;(0) + i(2)]_ K 2 [;(!);(!) _ i( 0)i(2)] J, 

1 



V 



(38) 
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FIG. 4. (Color online) Regular part of the interacting spin 
conductivity, a^.% 3 (uj), in units of (g/is) 2 /h within the lad- 
der approximation (see Fig. [2j for the spatially anisotropic 
but spin isotropic (A = 1) 2D Heisenberg model. Here 
uj = uj/Qmax where Q. ma x = 2,4(1 + rf)Sa{S)/h. The plot 
is for 77 = Jy/Jx values of 1 (black circles), 0.7 (red stars), 
and 0.5 (blue diamonds). The dashed lines are a guide to the 
eye. 
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1.0 



1.5 



2.0 



FIG. 5. (Color online) Regular part of the interacting spin 
conductivity, ^'(u), in units of (g/js) 2 /h within ladder ap- 
proximation for the spatial and spin anisotropic model. Here 
lu = u/Q. ma x where £l ma x = 2J X A(1 + r])Sa(S)/h. The plot 
is for 77 = 0.5 & A = 1 (blue diamonds), r\ = 0.5 & A = 1.001 
(red stars), 77 = 0.5 & A = 1.01 (green triangles), and 
77 = 1.5 & A = 1.01 (black circles). The dashed lines are 
a guide to the eye. 



zero temperature (since the zero frequency regular spin 
conductivity weight survives). 
B. Drude Weight: Anisotropy and Temperature 

In this section we investigate the finite temperature 
properties of the spin Drude weight, ~Ds{r], A,T). At fi- 
nite temperature AF long-range order cannot exist in the 
purely 2D Heiscnbcrg model. However it is known from 
a non linear sigma model analysis that the spin-spin cor- 
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relation length is rather long, upto to T~J/ks, where J 
is the exchange constant and is the Boltzmann fac- 
tor (37l - [39j . Therefore spin wave like excitations could 
exist in such a situation and our finite temperature anal- 
ysis would still be applicable. Furthermore, anisotropy 
and the presence of long-range exchange interactions in 
real material systems can help stabilize a magnetically or- 
dered AF state. Based on these assumptions we compute 
the spin Drude weight (within the collision less approxi- 
mation) which is a measure of the ballistic nature of the 
system. Using Eq. [231 we obtain the following expression 
for the spin Drude weight, DsiVi A, T) 

D ±= J x Sa(S) 2 ^ ( Mk^jk V 1_ 

n 2T(1 + 7 ? )A3 N 2f \ e k J sinh 2 §f, ' 

(40) 

For our model the spin Drude weight depends on three 
parameters - rj, A, and T. 

In Fig. [6] we show the dependence of the Drude 
weight, Ds(T, 77, A) on temperature, spatial, and spin 
anisotropy. The plot is for 77 = 1 & A = 1 (black circles), 
77 = 0.7 & A = 1 (red squares), 77 = 0.5 & A = 1 (blue 
up triangles), 77 = 1 & A = 1.001 (green rectangles), 
?7 = 1 & A = 1.01 (orange stars), and 77 = 1 & A = 1.1 
(brown down triangles). We find that with decreasing 77 
(increasing spatial anisotropy) the Drude weight is en- 
hanced at finite temperature. However increasing spin 
anisotropy suppresses the Drude weight at finite temper- 
atures. In the zero temperature limit, the Drude weight 
vanishes for both the spatial and spin anisotropic model. 
But finite temperatures encourage ballistic transport due 
to a non- vanishing value of the drude weight. Based on 
the above results we conclude that a spatially anisotropic 
2D quantum Heiscnbcrg AF is a better spin conductor 
than a spin anisotropic AF at finite temperatures. 

The conclusions for the non vanishing Drude weight 
at finite temperature should be taken with a caution- 
ary note. Within the collision less (that is non interact- 
ing) quasi particle approximation a finite Drude weight 
at non-zero temperatures is expected. In this approxi- 
mation the magnons propagate as undamped excitations 
with infinite lifetime and energy (40l . |4lj . However if 
collisions are included the magnon lifetime will become 
finite, the conductivity will no longer be ballistic, and 
will lead to a broadening of the delta function [4j], [42| ■ 
The broadening is determined by the inverse lifetime of 
the magnons. Although we do not perform the general 
analysis of the collision dominated transport [4(| to treat 
the singular part of the spin conductivity, we believe the 
qualitative effects of spatial and spin anisotropy on the 
Drude weight as predicted from the non-interacting ap- 
proximation will not be affected by the more rigorous 
collision dominated transport treatment. 

C. Collinear antiferromagnetic phase 

The spatially anisotropic Heiscnbcrg model can sup- 
port a CAF (stripe) phase - columnar or row. In this sec- 




T 



FIG. 6. (Color online) Temperature (T), spatial anisotropy 
parameter (77), and spin anisotropy parameter (A) depen- 
dence of the Drude weight, Dg(T, 77, A). The plot is for 
?7 = 1 & A = 1 (black circles), 77 = 0.7 & A = 1 (red squares), 
77 = 0.5 & A = 1 (blue up triangles), 77 = 1 & A = 1.001 
(green rectangles), 77 = 1 & A = 1.01 (orange stars), and 
?7 = 1& A = 1.1 (brown down triangles). The dashed lines 
are a guide to the eye. 



tion we explore the effects of spatial and spin anisotropy 
on the regular part of the spin conductivity in the pre- 
viously unexplored CAF phase. We solve the problem 
within the non-interacting approximation. The non- 
interacting limit itself presents interesting features in the 
spin conductivity (see Fig. [7]). To study the CAF phase 
we consider the following model Hamiltonian 

H = J x ^2 S?-Sf + -J y ^2 (Sf ■ Sf + Sf ■ Sf) 

<i,j> <i,j> 

i,y 

+ \{Sf + Sf +y + Sf + Sf + ~) + ASf 'O (41) 

The superscripts A and B refer to the up- and down- spin 
sublattices. The exchange interaction in the x-direction 
is denoted by J a and along the y-direction by J. y . The 
spin vector, S, is defined on every lattice site where the 
index i runs over all lattice sites and i+x and i +y runs 
over nearest neighbors in the x- and y- direction respec- 
tively. The S" 4,31 * 1 and S A,Bz operators have their usual 
meaning. The spin anisotropy parameter is given by A 
and the spatial anisotropy parameter by 77 = As in 
the AF phase the Hamiltonian can be rearranged as 

H = E + H U (42) 

where Eg is the classical energy and Hi is the quadratic 
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LSWT contribution (see [43|]) given by 



In the above we have the following 



1 - 



, «; k = A - ?y(A - cos k y ), (44) 
2J lx SekKk, 



C = A-r){A + l). 





0.7 


(43) 


0.6 




0.5 




3- 0.4 


(44) 


b 0.3 


(45) 


0.2 


(46) 


0.1 


(47) 


0.0 



Following the same procedure as in the AF phase we 
can define a spin current operator in the CAF phase and 
Fourier transform to obtain the DM representation 



Jxo 



2j\ x s 



s\n(k x ) B k 



£k«k ^4k 

K 

The coefficients Ak and Bk are given by 
ylk = A — ?/(A — cosky), £> k = 



(a+a k + /3+/3k). (48) 



COS k-r 



(49) 



We then obtain the expression for the regular part of the 
spin conductivity in the CAF phase as 



reg 

xx 



h 



C 2 N 



E 



jC 2 N 



E 



sin(k x ] 

£kK k 

sin(fc a ) 
£k«k 



1 



8{u 



2^). (50) 



where the symbols have been defined before. 

In Fig. [7] we display the results of the regular part 
of the interacting spin conductivity, <j xx 9 (uj), in units of 
(g/Lts) 2 /h, for the spatially anisotropic 2D Heiscnberg 
model in the CAF phase. The plot is for r\ = J y /J x 
values of -0.25 (red stars), -0.5 (blue diamonds), and -1 
(green triangles), -1.25 (orange rectangles), and -2 (black 
circles). The regular part does not vanish in the zero 
frequency limit for non-zero values of spatial anisotropy. 
Contrary to the AF phase in the CAF phase increasing 
spatial anisotropy leads to a decrease in the spin conduc- 
tivity. Also within the non-interacting approximation the 
spin conductivity in the CAF phase has a spike, a feature 
which is absent in the AF phase. The spin configuration 
in the CAF phase consists of ferromagnetic ordering in 
one direction and AF in the other. This is similar to a 
set of coupled 1D-AF chains in the case of a spatially 
anisotropic system. The sharp spike may be a signature 
of the one dimensional spin conductivity response fl9| . 

In Fig. [5] we show the results of the regular part of 
the non- interacting spin conductivity, a xx 9 (tu), in units 
of (g/^) 2 /h, for both the spatial and spin anisotropic 
2D Heiscnberg model in the CAF phase. The plot is for 
i) = Jy/Jx values of -0.25 k A=1.001 (red stars), -0.25 




FIG. 7. (Color online) Regular part of the non-interacting 
spin conductivity, a^% 3 (uj), in units of (g/^s) 2 /h, for the spa- 
tially anisotropic but spin isotropic 2D Heisenberg model 
in the collinear antiferromagnetic phase. The plot is for 
Tj — Jy/Jx values of -0.25 (red stars), -0.5 (blue diamonds), 
and -1 (green triangles), -1.25 (orange rectangles), and -2 
(black circles). The dashed lines are a guide to the eye. The 
regular part does not vanish in the spatially anisotropic model 
in the zero frequency limit for non-zero values of r\. 



k A=1.01 (black triangles), and -0.25 k A=l.l (blue 
circles). As in the AF phase, inclusion of the smallest 
amount of spin anisotropy suppresses the zero frequency 
value of <j xx 9 {uj). The gap in the spin conductivity is 
controlled by spin anisotropy only. Finally, similar to 
the AF phase, we find that the spatially anisotropic CAF 
phase magnets are better spin conductors than those with 
spin anisotropy. 




FIG. 8. (Color online) Regular part of the non-interacting 
spin conductivity, a^% s (u)), in units of (g/^s) 2 /h, for the spa- 
tial and spin isotropic 2D Heisenberg model in the collinear 
antiferromagnetic phase. The plot is for 77 = J y /J x values of 
r? = -0.25 & A = 1.001 (red stars), r) = -0.25 & A = 1.01 
(black triangles), and 77 = —0.25 & A = 1.1 (blue circles). 
The dashed lines are a guide to the eye. Inclusion of the 
smallest amount of spin anisotropy A introduces a gap. This 
is similar to the spin conductivity behavior in the AF phase. 
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V. SUMMARY AND CONCLUSION 

In this paper we analyze the effect of spatial and spin 
anisotropy on spin conductivity in a 2D S=l/2 Hciscn- 
berg quantum magnet on a square lattice. We explore 
the model in both the Neel AF phase and the CAF 
phase. Based on a many body Green function formal- 
ism and utilizing the Kubo formula for spin transport 
we explore the effects of spin and spatial anisotropy in 
the Hcisenberg model. We utilize material parameters 
from SrZnVO(P04)2 (spatially anisotropic) and La2Ni04 
(spin isotropic). We find that these anisotropics have 
opposite effects on the magnetic system. In the AF 
phase spatial anisotropy does not suppress the regular 
part of the spin conductivity in the zero frequency limit 
- rather it enhances it. However spin anisotropy re- 
duces the weight. In the CAF phase (within the non- 
interacting approximation) spin conductivity has a fi- 
nite value for non-zero values of spatial anisotropy pa- 
rameter which decreases as the spatial anisotropy is en- 
hanced. Furthermore, we find that the CAF phase dis- 
plays a spike in the spin conductivity not seen in the AF 
phase. The spike could be a signature of the quasi-lD 
behavior of the spatially anisotropic system. In both the 
AF and CAF phase we find that inclusion of the small- 
est amount of spin anisotropy causes o- reg (uj) to develop 
a gap. We also explore the finite temperature Drudc 
weight, Ds(t?, A, T), in the AF phase for various spatial 
and spin anisotropy parameters. We find that within the 
collision less (non-interacting) quasi particle approxima- 
tion the Drude weight goes to zero as the temperature 
approaches zero and remains non-zero at finite tempera- 
tures. Increasing spatial anisotropy increases the Drude 
weight value and increasing spin anisotropy decreases the 



Drudc weight value. Based on these studies we conclude 
that quantum Heisenberg antiferromagnets with spatial 
anisotropy are better spin conductors than those with 
spin anisotropy at both zero and finite temperatures. 
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Appendix A: DM Interaction Vertices 

In this appendix wc list the quartic interaction vertices 
obtained in the AF phase after the DM transformation of 
the orginal Hamiltonian, Eq. [TJ The symbols 1, 2, 3, and 
4 stand for wavevectors ki, k2, k3, and k4 respectively. 
See Eq. [9] for the definition of Uk, Vk, and Eq[TT]for £k- 
In the vertices below wc use the following definitions 



[^1234 = W1U2W3U4, 

1/2 



ik = 



"k 



1 



(Al) 
(A2) 



U (1) = -£/i234{ii [14(137134 - A714) - (AX3713 - 71)] + x 2 [x 4 (1-37234 - A724) - (A13723 - 72)]}, (A3) 

U (2) = -2[/i234{27£2[£4(A7i4 - 137134) - (ji - Ax 3 7i 3 )] + [x 4 (A7 2 4 - ^37234) - (72 - Al37 23 )]}, ( A4 ) 

U (3) = -2[/i23 4 {2;i[(A7i3 - 1371) ~ 14(7134 - A13714)] + i 2 [(A723 - 1372) - 14(7234 - Ai 3 7 24 )]}, (A5) 

U (4) = -4t r i 2 3 4 {lll 2 [(l37l34 - A714) - l 4 (Al 3 7i3 - 71)] + [(l 3 7234 ~ A7 24 ) - l 4 (Al 3 7 2 3 - 72)]}, (A6) 

U (5) = -2t r i 234 {i2[i4(Ai37i3 - 71 ) - (137134 - A714)] + xi[x 4 (Ai3723 - 72) - (137234 - A7 24 )]}, (A7) 

U (6) = -2£/i 2 34{!i! 2 [(A!3724 " 7234) ~ ^4(^372 ~ A723)] + [(Al 2 7i 4 - 7134 ) ~ l4(l37l ~ A713)]}, (A8) 

U (7) = -£/i23 4 {ii[(7i34 - A13714) - i 4 (A7i 3 - 1371)] + 12K7234 - A13724) - x 4 (A72 3 - 1372)]}, (A9) 

U (8) = -J7i 2 34{l 2 [l4(l37l34 ~ A714) - (A13713 - 71)] + li[l 4 (l37234 - A724) - (A13723) - 72]}, (A10) 

U (9) = -J7i 2 34{i 2 [(7i34 - A13714) - x 4 (A7i3 - 1371)] + ii[(7234 - A13724) - (i 4 (A7 23 - 1372)]}- (All) 
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